Dichotomy between large local and small ordered magnetic moment in Iron-based 

superconductors 
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We study a four band model for iron-based superconductors within local density approximation 
+ dynamical mean field theory (LDA+DMFT). This successfully reproduces the results of models 
which take As p degrees of freedom explicitly into account and has several physical advantages 
over the standard five d-band model. Our findings reveal that the new superconductors are more 
strongly correlated than their single-particle properties suggest. Two-particle correlation functions 
unveil the dichotomy between local and ordered magnetic moments in these systems, calling for 
further experiments to better resolve the short time scale spin dynamics. 
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In the recently discovered iron-based superconductors 
[l[ the role of electronic correlation is still highly un- 
clear. Strongly correlated materials are characterized by 
the presence of large local magnetic moments, which typi- 
cally long-range order if temperatures are sufficiently low. 
In the proximity of such magnetic phases, also supercon- 
ductivity is often observed, with the close-by magnetic 
fluctuations usually cited as evidence for unconventional 
(not phonon-mediated) superconductivity. If the local 
magnetic moment is small and the system is metallic, 
weak-coupling theories like local spin density approxi- 
mation (LSDA) can be applied. On the other hand, 
if the local magnetic moment is large and the exchange 
coupling between neighboring spins is the dominant in- 
teraction, not only the electronic states around the Fermi 
level but also the higher energy excitations (on the scale 
of the local Hubbard interaction U) are expected to play 
a role in the superconducting pairing mechanism. For in- 
stance, the latter is definitely the case for cuprates, which 
are Mott insulating in the absence of carrier doping. Iron 
pnictides instead are metallic and undergo a spin-density 
wave (SDW) transition below r?a 150K whose character- 
istics are still under debate. Understanding the nature of 
the magnetic properties can therefore also help to clarify 
the origin of superconductivity in these materials. 

Experimentally it has been clarified that the differ- 
ent members of the pnictide family have quite differ- 
ent ordered magnetic moments Q, ranging from O.S/is 
(or d) in LaFeAsO (1111 compound) to 2.2^b 

in FeTe (11 compound). The band structures of these 
compounds however do not show distinctive differences, 
and indeed, LSDA always yields an ordered moment of 
~2.0^B, for the experimental crystal structures 
The failure of density functional theory (DFT) in cap- 
turing the correct ground-state properties of pnictides is 
hence a very important point which needs to be carefully 
analyzed in order to understand the physics of these com- 



pounds: While usually LSDA underestimates the size of 
the ordered moment, here the opposite happens. Even 
more important than the size of the ordered moment is 
the fact that in these systems the magnetic properties 
are extremely sensitive to the choice of the exchange- 
correlation functional or of the crystal structure In 
this situation, and also with the small magnetic moment 
indicating the proximity to a quantum-critical point, 
quantum fiuctuations strongly influence the physics and 
the moments of iron-based superconductors. It is there- 
fore natural to conclude that dynamical quantum fluc- 
tuations, not included in LSDA, are crucial for these 
systems. In particular they can explain the presence of 
large local magnetic moments which form because of local 
Coulomb and exchange interaction but only give rise to 
a much smaller ordered moment at lower temperatures. 

There have been many attempts to go beyond LSDA 
taking electronic correlations more accurately into ac- 
count. Among these, dynamical mean field theory 
(DMFT) is one of the most promising, particularly when 
combined with ab-initio band structure calculations 
However, the results of such LDA-I-DMFT calculations 
for iron-based superconductors @- 13| strongly depend on 



which orbitals are included in the ab-initio one-particle 
Hamiltonian and on the values of the interaction param- 
eters used. As a consequence, DMFT calculations have 
been employed by different groups in fairly different ways, 
namely to support that iron pnictides are strongly, inter- 
mediately, or weakly correlated, respectively. The ma- 
jority of these studies focused on single-particle spectra 
and on the comparison with photoemission experiments, 
except for Refs. |lj and ll3| where also the spin suscepti- 
bility has been calculated. In this paper we focus on the 
dynamics of the local magnetic moment, which we argue 
is a key indicator for understanding the physics of iron- 
based superconductors. In particular, we conclude that 
in the single-particle spectral function correlation effects 
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are hardly visible, while, at the same time, the spin-spin 
correlation function reveals the existence of a large local 
magnetic moment. This turns out to be crucial for the 
explanation of some controversial experimental results in 
these systems. 

For iron-based superconductors, two classes of mod- 
els have been proposed 15, 3]: One is a d only model 
which takes the Fe 3d degrees of freedom into account, 
while the others are dp or dpp models considering pnic- 
togen/chalcogen p and O 2p electrons explicitly. In this 
work, we take the d model as a starting point. The d 
models considered hitherto, however, pose some physi- 
cal and technical problems: Each Wannier function of d 
character has a fairly different spread in real space, due 
to the orbital-dependent hybridization between p and d 
15| . Thus the interaction parameters strongly depend on 
the orbital [l7, 18 1. In our coordinate system, in which x- 
and y-axes point to the pnictogen/chalcogen atom, the 
Sz^-r^ orbital, for example, has a small p-d hybridiza- 
tion, so that it is well localized. On the other hand, the 
x^-y^ orbital has long tails in the direction of the pnicto- 
gen or chalcogen sites. Such an orbital dependence causes 
problems when calculating the self-energy due to electron 
correlation: In order to avoid a double-counting of corre- 
lation effects already considered within LDA, one would 
have to introduce an ad hoc orbital-dependent level shift 
in the many-body calculation. This level shift has been 
realized to be particularly important for the Sz^-r^ or- 
bital. In calculations based on the so-called FLEX ap- 
proximation [l^ , it has been shown that the 3z^-r^ level 
becomes higher in energy and makes a large Fermi sur- 
face not present in LDA. A similar tendency is also seen 
in DMFT calculations , namely the 3z^-r^ occupancy 
gets dramatically smaller than in LDA depending on the 
strength of interaction parameters, contrary to what hap- 
pens in LDA-I-DMFT calculations for dp and dpp models 
[10l4l2| . In order to overcome this problem, some au- 
thors added a constant part to the self-energy [l^, or 
constrained the zero frequency value of the self-energy to 
get the appropriate orbital shift [i^l . 



Here we take another route and assume the 3z 



to 



be fully occupied. Wc therefore do not include it in 
our low-energy Hamiltonian for LaFcAsO and arc left 
with a model with merely four equally correlated orbitals 
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which allows us to circumvent an orbital dependent 
double-counting adjustment. This approximation is jus- 
tified by the fact that the band with mainly Sz^-r^ char- 
acter lies below the Fermi level. We argue that the results 
of the four orbital model compare to the dpp model much 
better than the five band one. It has been already shown 
2lj that the Fermi surface of the four band model is al- 
most exactly the same as that of the five band model in 
LDA, and docs not change so drastically even after the 
inclusion of many-body effects. 

The four and the five d-band models have of course dif- 
ferent values of the interaction parameters. Constrained 
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FIG. 1: Orbital-resolved spectral function of the four-band 
model for LaFeAsO at T = 460K, for U = 1.8eV, J = 0.45eV 
(solid lines) and U = J = (dashed lines). Including electronic 
correlations does not change the spectrum drastically. 



random phase approximation (RPA) calculations for the 
five band model give an intra-orbital Coulomb interac- 
tion U of about 2.2 — 3.3eV and a Hund's rule coupling 
of about 0.3 — 0.6eV . In our case we have to consider 
the screening effect of 3z^-r^ orbital as well as the slightly 
more extended tails of the Wannier functions which both 
reduce the value of U. For this reason, hereafter we use 
U = 1.8eV. Numerical limitations necessitate a Hund's 
exchange J of Ising-type for four orbitals. Below, we will 
show however a comparison of results between Ising and 
full SU(2) symmetric interaction, for a model with fewer 
orbitals. The differences turn out not to be relevant for 
the present discussion. 

In Fig. [TJ we show the spectral functions of the four 
band model for non-interacting electrons (dashed lines) 
and for C/ = 1.8eV, J = 0.45eV and temperature T = 460K 
(/3 = 25cV^^). Including U and J renormalizes (shrinks 
the width of) the minimum structure around the Fermi 
level. At the same time, some of the spectral weight is 
shifted to Hubbard-like shoulders at higher energies, in 
a way such that the overall bandwidth and spectrum re- 
mains close to the non-interacting one. The values of 
the quasiparticle weights are Z ~0.51, 0.45, and 0.60 for 
yz(xz), x^-y^, and xy band, respectivel y, in ag reement 
with photoemission spectroscopy (PES) [22|, |23|. As an- 
ticipated above, the results of the four orbital model re- 
sembles very closely the dp and the dpp models [ll|, [l3 ■ 
For instance, in our xy band the peak around '~ leV 



3 



shifts towards lower energies, very similarly to the corre- 
sponding spectral function in Ref. 12 and Ref. [Tl| (there 



denoted as x^-y'^ orbital as their coordinate system is 
rotated by 45°). On the other hand, our x^-y'^ has a 
structure around ~ leV coming form the hybridization 



to the 2>z^-r^ and the pnictogen/chalcogen p, which are ^ 



included effectively in our four band Hamiltonian. We 
can therefore conclude that the four band model reliably 
reproduces the results of photoemission experiments with 
the clear advantage of (i) having a smaller number of pa- 
rameters and (ii) of yielding a set of d orbitals with much 
more similar spatial spread (and double counting correc- 
tion). 

The picture arising from merely analyzing the spectral 
function hence suggests that iron-based superconductors 
are quite far from being standard strongly correlated ma- 
terials, such as cuprates or other transition-metal com- 
pounds. On the other hand, calculations based on FLEX 
which one would expect to work for weakly correlated 
materials, here fail to reproduce the correct stripe pat- 
tern of antiferromagnetic spin fluctuations (2ll |. This in- 
dicates that iron-based superconductors cannot be cate- 
gorized as weakly correlated systems, at least not with 
respect to their two-particle correlation functions. 

These considerations naturally leads to the question: 
Are iron-based superconductors more correlated than 
their single particle quantities such as the photoemis- 
sion spectra suggest? To answer this question we cal- 
culated the local spin susceptibility of LaFeAsO within 
LDA-I-DMFT and study whether or not this indicates 
the existence of a large local magnetic moment in these 
compounds. 

In Fig. [21 wc plot the (dynamical) local spin-spin cor- 
relation function Xim (r) = (5f (t)5^,(0)) for (imaginary) 
time T. Resolved are its intra-orbital (I = m) and inter- 
orbital {^i^jnXim{T)) contribution. Similar to the case 
of one-particle properties, the [/-driven intra-orbital spin 
correlation is only slightly enhanced in comparison to 
the non-interacting value which for equal times (r = 0) 
is [23l- In stark contrast, the inter-orbital con- 

tribution which vanishes without interaction is strongly 
enhanced. This reflects the strong tendency of the sys- 
tem to align spins between different orbitals. It can be 
understood by noting that, since the crystal field split- 
ting is small (^0.2cV) [2l|, even an intermediate value of 
the Hund's rule exchange J is very effective. Hence, the 
large inter-orbital spin correlation function points to the 
importance J plays for inducing electronic correlations 
in iron-based superconductors. 

Let us note in passing the inset of Fig. [5] which com- 
pares the total xi^) for a half-filled two-band Hubbard 
model with a semi-elliptic density of states with band- 
width 4eV, U = 2.5eV, J = 0.5eV and /? = 20eV-\ 
a choice that gives renormalization factors very similar 
to our four-band realistic calculation. As already men- 
tioned, the close agreement of the two curves allows us 
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FIG. 2: Spin-spin correlation function for J = 0.45eV and 
/3 = 25eV~^. We plot the different intra-orbital contributions 
and the sum of all inter-orbital contribution. This orbital- 
resolved presentation clearly shows that the J-induced inter- 
orbital correlation is particularly large. Inset: Comparison 
between Ising- and SU(2)-symmetric Hund's exchange for a 
related two-band model, showing only quite small differences. 
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FIG. 3: Total spin-spin correlation function for LaFeAsO at 
two different values of J and /3 = 25eV~^, compared to the 
non-interacting U = J = case. The short-time (r — 0) 
LDA+DMFT local moment for J = 0.45eV is mioc = 2.16^is , 
comparably large as in LSDA; whereas the long-time moment 
is screened to only m^Q.l^B at r = 50K, see text. 



to exclude that the Ising approximation for the J term 
has a big influence in the relevant parameter regime. 

Turning to the central result of our paper, the total 
x(t) in Fig. [3l we see a much larger susceptibility than 
in the non-interacting case. This enhancement of x(t) is 
a measure for the correlations in iron-based superconduc- 
tors. At T = 0, the susceptibihty x(r) gives us a direct 
measure of the bare local moment toioc = ■\/x(0). This 
local moment corresponds to the responses on short time 
(or high energy) scales and is quite large, i.e., mioc = 
2.16^B and 2.45/i_B for J = 0.45eV and 0.5 eV, respec- 
tively. Note that such variations of J are realistic when 



4 



going from LaFeAsO to FeTe. While the J = 0.45 eV 
value of TOioc is only slightly T-dependent (TOioc = 2.36/iB 
at /? = 10 eV^^), at J = 0.5 eV an even larger mo- 
ment is formed at higher temperature (mioc = 3.44/xb 
at /3 = 10eV-i). 

While this local moment is large and indicates strong 
correlations, it is not the one which was hitherto mea- 
sured experimentally. Experiments such as magnetic sus- 
ceptibility measurements, nuclear magnetic resonance, 
Mofibauer spectroscopy and muon relaxation are slow 
compared to the electronic dynamics on the fs time scale. 
Hence, these experiments correspond to larger r's or the 
integrated (static) susceptibility x(w = 0) = J^^ drxir). 
A central result of our calculation is that this long- 
time susceptibility or a corresponding magnetic moment, 
which one can define by resolving x(uj = 0) = m? jT 
for m, is strongly reduced (screened) compared to the 
instantaneous local magnetic moment m\oc- Already at 
/3 = 25eV~^ (/3 = 10eV~^) the dynamic screening leads 
to strongly reduced moments of m = 1.2(1.9) and 1.8(3.4) 
for J ~ 0.45 eV and 0.5 eV, respectively. And for lower 
temperatures these values are much further reduced be- 
cause of the screening. At T = 50 K, i.e., in the temper- 
ature range where such magnetic moments were experi- 
mentally measured, an extrapolation of our data yields 
a crude (overestimated) approximation of m « 0.7^b 
(25I [26j . Since also the SDW phase of iron-based su- 
perconductors is very itinerant for both spin species, we 
expect that the dynamical screening identified here as 
the origin of the smallness of the (long-time) magnetic 
moment survive (to a large extent) also in the magnetic 
phase. 

In conclusion, LDA-I-DMFT predicts that the local 
magnetic moment in iron-based superconductor is, in 
the paramagnetic phase, comparable to the ordered mo- 
ment of LSDA. This moment is formed due to a local 
Hund's rule spin alignment. However, there is a di- 
chotomy between this local magnetic moment and the 
dynamically screened moment, which is much smaller 
and beyond LSDA. The latter is actually much smaller 
in LDA-I-DMFT. Experiments performed hitherto mea- 
sured the low-energy (or long-time) moment, i.e., the dy- 
namically screened one. For measuring the (bare) local 
moment, experimental measurements on the time scale 
of fs are needed. A possibility to this end is to integrate 
neutron scattering measurements over Q and lu. For such 
an experiment, our calculations predict an intermediate- 
to-large value of the local magnetic moment. Similarly, 
X-ray spectroscopy is a very promising technique for mea- 
suring the size of the local magnetic moment, but so far 
this has been mainly used to estimate the strength of 
the interaction parameter, like, e.g., in Ref. |271 Such 
experiments, if performed, will clarify whether our idea 
of iron-based superconductors being more strongly corre- 
lated than what is naively expected from photocmission 
experiments, is correct. This can eventually settle the 



role electronic correlations play in the new class of iron- 
based superconductors. 
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